Association studies between COVID-19 and SSc-ILD

Severe COVID-19 patients may develop pulmonary fibrosis, similar to SSc-ILD disease, suggesting a potential link between the two diseases. However, there are limited treatment options for SSc-ILD-type diseases. Therefore, investigating pathological markers of the two diseases can provide valuable insights for treating related conditions. RNA sequencing technology offers high throughput and precision. However, the bimodal nature of RNA-Seq data cannot be accurately captured by commonly used algorithms such as DESeq2. To address this issue, the Beta-Poisson model has been developed to identify differentially expressed genes. Unlike the classical DESeq2 algorithm, the Beta-Poisson model introduces a Beta distribution to construct a new hybrid distribution in place of the Gamma distribution of the Gamma-Poisson distribution, effectively characterizing the bimodal features of RNA-Seq data. The transcriptomes of SARS-CoV infection and SSc-ILD disease in the lung epithelial cell dataset were analyzed to identify common differentially expressed genes of SARS-CoV and SSc-ILD disease. Gene function and signaling pathway enrichment analysis and protein-protein interaction (PPI) network were used to identify common pathways and drug targets for SSc-ILD with COVID-19 infection. The results show that there are 50 differentially expressed genes in common between COVID-19 and SSC-ILD. The functions of these genes are mainly enriched in immune system response, interferon signaling pathway and other related signaling pathways, and enriched in biological processes such as cell defense response to virus and interferon regulation. Based on the detection of hub genes based on PPIs network, it is predicted that STAT1, ISG15, IRF7, MX1, EIF2AK2, DDX58, OAS1, OAS2, IFIT1 and IFIT3 are the key genes involved in the pathological phenotype of the two diseases. Based on the key genes, the interaction of transcription factor (TF) and miRNA with common differentially expressed genes is also identified. The possible pathological markers of the two diseases and related molecular regulatory mechanisms of disease treatment are revealed to provide theoretical basis for the treatment of the two diseases. © 2023 Editorial Office of Journal of Shenzhen University. All rights reserved.

An efficient computational protocol for template-based design of peptides that inhibit interactions involving SARS-CoV-2 proteins.

The RNA-dependent RNA polymerase (RdRp) complex of SARS-CoV-2 lies at the core of its replication and transcription processes. The interfaces between holo-RdRp subunits are highly conserved, facilitating the design of inhibitors with high affinity for the interaction interface hotspots. We, therefore, take this as a model protein complex for the application of a structural bioinformatics protocol to design peptides that inhibit RdRp complexation by preferential binding at the interface of its core subunit nonstructural protein, nsp12, with accessory factor nsp7. Here, the interaction hotspots of the nsp7-nsp12 subunit of RdRp, determined from a long molecular dynamics trajectory, are used as a template. A large library of peptide sequences constructed from multiple hotspot motifs of nsp12 is screened in-silico to determine sequences with high geometric complementarity and interaction specificity for the binding interface of nsp7 (target) in the complex. Two lead designed peptides are extensively characterized using orthogonal bioanalytical methods to determine their suitability for inhibition of RdRp complexation. Binding affinity of these peptides to accessory factor nsp7, determined using a surface plasmon resonance (SPR) assay, is slightly better than that of nsp12: dissociation constant of 133nM and 167nM, respectively, compared to 473nM for nsp12. A competitive ELISA is used to quantify inhibition of nsp7-nsp12 complexation, with one of the lead peptides giving an IC50 of 25µM . Cell penetrability and cytotoxicity are characterized using a cargo delivery assay and MTT cytotoxicity assay, respectively. Overall, this work presents a proof-of-concept of an approach for rational discovery of peptide inhibitors of SARS-CoV-2 protein-protein interactions.

Structure and Functions of HMGB2 Protein.

High-Mobility Group (HMG) chromosomal proteins are the most numerous nuclear non-histone proteins. HMGB domain proteins are the most abundant and well-studied HMG proteins. They are involved in variety of biological processes. HMGB1 and HMGB2 were the first members of HMGB-family to be discovered and are found in all studied eukaryotes. Despite the high degree of homology, HMGB1 and HMGB2 proteins differ from each other both in structure and functions. In contrast to HMGB2, there is a large pool of works devoted to the HMGB1 protein whose structure-function properties have been described in detail in our previous review in 2020. In this review, we attempted to bring together diverse data about the structure and functions of the HMGB2 protein. The review also describes post-translational modifications of the HMGB2 protein and its role in the development of a number of diseases. Particular attention is paid to its interaction with various targets, including DNA and protein partners. The influence of the level of HMGB2 expression on various processes associated with cell differentiation and aging and its ability to mediate the differentiation of embryonic and adult stem cells are also discussed.

Proteomic understanding of SARS-CoV-2 infection and COVID-19: Biological, diagnostic, and therapeutic perspectives

The research community responded rapidly to the SARS-CoV-2 pandemic with a burst of proteomic studies to understand this new virus. We focus in this chapter on proteomic approaches and advances in the areas of the viral proteome and the viral-host protein-protein interactions exploited to facilitate the pathogenic life cycle of SARS-CoV2. We also outline the proteome and posttranslational protein modifications of infected cells and blood proteomics of COVID-19 patients in relation to pathogenesis, host response, and disease severity. Finally, tools and application for COVID-19 proteomics and implementations in diagnostics and therapies are surveyed before summarizing the present achievements and future perspectives. © 2023 Elsevier Inc. All rights reserved.

Myxobacterial depsipeptide chondramides interrupt SARS-CoV-2 entry by targeting its broad, cell tropic spike protein.

The severity of the COVID-19 pandemic has necessitated the search for drugs against SARS-CoV-2. In this study, we explored via in silico approaches myxobacterial secondary metabolites against various receptor-binding regions of SARS-CoV-2 spike which are responsible in recognition and attachment to host cell receptors mechanisms, namely ACE2, GRP78, and NRP1. In general, cyclic depsipeptide chondramides conferred high affinities toward the spike RBD, showing strong binding to the known viral hot spots Arg403, Gln493 and Gln498 and better selectivity compared to most host cell receptors studied. Among them, chondramide C3 (1) exhibited a binding energy which remained relatively constant when docked against most of the spike variants. Chondramide C (2) on the other hand exhibited strong affinity against spike variants identified in the United Kingdom (N501Y), South Africa (N501Y, E484K, K417N) and Brazil (N501Y, E484K, K417T). Chondramide C6 (9) showed highest BE towards GRP78 RBD. Molecular dynamics simulations were also performed for chondramides 1 and 2 against SARS-CoV-2 spike RBD of the Wuhan wild-type and the South African variant, respectively, where resulting complexes demonstrated dynamic stability within a 120-ns simulation time. Protein-protein binding experiments using HADDOCK illustrated weaker binding affinity for complexed chondramide ligands in the RBD against the studied host cell receptors. The chondramide derivatives in general possessed favorable pharmacokinetic properties, highlighting their potential as prototypic anti-COVID-19 drugs limiting viral attachment and possibly minimizing viral infection.Communicated by Ramaswamy H. Sarma.

SARS-CoV-2 SUD2 and Nsp5 Conspire to Boost Apoptosis of Respiratory Epithelial Cells via an Augmented Interaction with the G-Quadruplex of BclII.

The molecular mechanisms underlying how SUD2 recruits other proteins of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to exert its G-quadruplex (G4)-dependent pathogenic function is unknown. Herein, Nsp5 was singled out as a binding partner of the SUD2-N+M domains (SUD2core) with high affinity, through the surface located crossing these two domains. Biochemical and fluorescent assays demonstrated that this complex also formed in the nucleus of living host cells. Moreover, the SUD2core-Nsp5 complex displayed significantly enhanced selective binding affinity for the G4 structure in the BclII promoter than did SUD2core alone. This increased stability exhibited by the tertiary complex was rationalized by AlphaFold2 and molecular dynamics analysis. In line with these molecular interactions, downregulation of BclII and subsequent augmented apoptosis of respiratory cells were both observed. These results provide novel information and a new avenue to explore therapeutic strategies targeting SARS-CoV-2. IMPORTANCE SUD2, a unique protein domain closely related to the pathogenesis of SARS-CoV-2, has been reported to bind with the G-quadruplex (G4), a special noncanonical DNA structure endowed with important functions in regulating gene expression. However, the interacting partner of SUD2, among other SARS-CoV-2 Nsps, and the resulting functional consequences remain unknown. Here, a stable complex formed between SUD2 and Nsp5 was fully characterized both in vitro and in host cells. Moreover, this complex had a significantly enhanced binding affinity specifically targeting the Bcl2G4 in the promoter region of the antiapoptotic gene BclII, compared with SUD2 alone. In respiratory epithelial cells, the SUD2-Nsp5 complex promoted BclII-mediated apoptosis in a G4-dependent manner. These results reveal fresh information about matched multicomponent interactions, which can be parlayed to develop new therapeutics for future relevant viral disease.

Modeling of SARS-CoV-2 Virus Proteins: Implications on Its Proteome.

COronaVIrus Disease 19 (COVID-19) is a severe acute respiratory syndrome (SARS) caused by a group of beta coronaviruses, SARS-CoV-2. The SARS-CoV-2 virus is similar to previous SARS- and MERS-causing strains and has infected nearly six hundred and fifty million people all over the globe, while the death toll has crossed the six million mark (as of December, 2022). In this chapter, we look at how computational modeling approaches of the viral proteins could help us understand the various processes in the viral life cycle inside the host, an understanding of which might provide key insights in mitigating this and future threats. This understanding helps us identify key targets for the purpose of drug discovery and vaccine development.

Pathogen-driven cancers from a structural perspective: Targeting host-pathogen protein-protein interactions.

Host-pathogen interactions (HPIs) affect and involve multiple mechanisms in both the pathogen and the host. Pathogen interactions disrupt homeostasis in host cells, with their toxins interfering with host mechanisms, resulting in infections, diseases, and disorders, extending from AIDS and COVID-19, to cancer. Studies of the three-dimensional (3D) structures of host-pathogen complexes aim to understand how pathogens interact with their hosts. They also aim to contribute to the development of rational therapeutics, as well as preventive measures. However, structural studies are fraught with challenges toward these aims. This review describes the state-of-the-art in protein-protein interactions (PPIs) between the host and pathogens from the structural standpoint. It discusses computational aspects of predicting these PPIs, including machine learning (ML) and artificial intelligence (AI)-driven, and overviews available computational methods and their challenges. It concludes with examples of how theoretical computational approaches can result in a therapeutic agent with a potential of being used in the clinics, as well as future directions.

SARS-COV-2 Coronavirus Papain-like Protease PLpro as an Antiviral Target for Inhibitors of Active Site and Protein-Protein Interactions.

The papain-like protease PLpro of the SARS-CoV-2 coronavirus is a multifunctional enzyme that catalyzes the proteolytic processing of two viral polyproteins, pp1a and pp1ab. PLpro also cleaves peptide bonds between host cell proteins and ubiquitin (or ubiquitin-like proteins), which is associated with a violation of immune processes. Nine structures of the most effective inhibitors of the PLpro active center were prioritized according to the parameters of biochemical (IC 50) and cellular tests to assess the suppression of viral replication (EC 50) and cytotoxicity (CC 50). A literature search has shown that PLpro can interact with at least 60 potential protein partners in cells, 23 of which are targets for other viral proteins (human papillomavirus and Epstein-Barr virus). The analysis of protein-protein interactions showed that the proteins USP3, UBE2J1, RCHY1, and FAF2 involved in deubiquitinylation and ubiquitinylation processes contain the largest number of bonds with other proteins; the interaction of viral proteins with them can affect the architecture of the entire network of protein-protein interactions. Using the example of a spatial model of the PLpro/ubiquitin complex and a set of 154 naturally occurring compounds with known antiviral activity, 13 compounds (molecular masses in the range of 454-954 Da) were predicted as potential PLpro inhibitors. These compounds bind to the "hot" amino acid residues of the protease at the positions Gly163, Asp164, Arg166, Glu167, and Tyr264 involved in the interaction with ubiquitin. Thus, pharmacological effects on peripheral PLpro sites, which play important roles in binding protein substrates, may be an additional target-oriented antiviral strategy.

Investigating the competition between ACE2 natural molecular interactors and SARS-CoV-2 candidate inhibitors.

The SARS-CoV-2 pandemic still poses a threat to the global health as the virus continues spreading in most countries. Therefore, the identification of molecules capable of inhibiting the binding between the ACE2 receptor and the SARS-CoV-2 spike protein is of paramount importance. Recently, two DNA aptamers were designed with the aim to inhibit the interaction between the ACE2 receptor and the spike protein of SARS-CoV-2. Indeed, the two molecules interact with the ACE2 receptor in the region around the K353 residue, preventing its binding of the spike protein. If on the one hand this inhibition process hinders the entry of the virus into the host cell, it could lead to a series of side effects, both in physiological and pathological conditions, preventing the correct functioning of the ACE2 receptor. Here, we discuss through a computational study the possible effect of these two very promising DNA aptamers, investigating all possible interactions between ACE2 and its experimentally known molecular partners. Our in silico predictions show that some of the 10 known molecular partners of ACE2 could interact, physiologically or pathologically, in a region adjacent to the K353 residue. Thus, the curative action of the proposed DNA aptamers could recruit ACE2 from its biological functions.

Timeline of changes in spike conformational dynamics in emergent SARS-CoV-2 variants reveal progressive stabilization of trimer stalk with altered NTD dynamics.

SARS-CoV-2 emergent variants are characterized by increased viral fitness and each shows multiple mutations predominantly localized to the spike (S) protein. Here, amide hydrogen/deuterium exchange mass spectrometry has been applied to track changes in S dynamics from multiple SARS-CoV-2 variants. Our results highlight large differences across variants at two loci with impacts on S dynamics and stability. A significant enhancement in stabilization first occurred with the emergence of D614G S followed by smaller, progressive stabilization in subsequent variants. Stabilization preceded altered dynamics in the N-terminal domain, wherein Omicron BA.1 S showed the largest magnitude increases relative to other preceding variants. Changes in stabilization and dynamics resulting from S mutations detail the evolutionary trajectory of S in emerging variants. These carry major implications for SARS-CoV-2 viral fitness and offer new insights into variant-specific therapeutic development.

Liaisons dangereuses: Intrinsic Disorder in Cellular Proteins Recruited to Viral Infection-Related Biocondensates.

Liquid-liquid phase separation (LLPS) is responsible for the formation of so-called membrane-less organelles (MLOs) that are essential for the spatio-temporal organization of the cell. Intrinsically disordered proteins (IDPs) or regions (IDRs), either alone or in conjunction with nucleic acids, are involved in the formation of these intracellular condensates. Notably, viruses exploit LLPS at their own benefit to form viral replication compartments. Beyond giving rise to biomolecular condensates, viral proteins are also known to partition into cellular MLOs, thus raising the question as to whether these cellular phase-separating proteins are drivers of LLPS or behave as clients/regulators. Here, we focus on a set of eukaryotic proteins that are either sequestered in viral factories or colocalize with viral proteins within cellular MLOs, with the primary goal of gathering organized, predicted, and experimental information on these proteins, which constitute promising targets for innovative antiviral strategies. Using various computational approaches, we thoroughly investigated their disorder content and inherent propensity to undergo LLPS, along with their biological functions and interactivity networks. Results show that these proteins are on average, though to varying degrees, enriched in disorder, with their propensity for phase separation being correlated, as expected, with their disorder content. A trend, which awaits further validation, tends to emerge whereby the most disordered proteins serve as drivers, while more ordered cellular proteins tend instead to be clients of viral factories. In light of their high disorder content and their annotated LLPS behavior, most proteins in our data set are drivers or co-drivers of molecular condensation, foreshadowing a key role of these cellular proteins in the scaffolding of viral infection-related MLOs.

SARS-CoV-2 ORF8 dimerization and binding mode analysis with class I MHC: computational approaches to identify COVID-19 inhibitors.

SARS-CoV-2 encodes eight accessory proteins, one of which, ORF8, has a poorly conserved sequence with SARS-CoV and its role in viral pathogenicity has recently been identified. ORF8 in SARS-CoV-2 has a unique functional feature that allows it to form a dimer structure linked by a disulfide bridge between Cys20 and Cys20 (S-S). This study provides structural characterization of natural mutant variants as well as the identification of potential drug candidates capable of binding directly to the interchain disulfide bridge. The lead compounds reported in this work have a tendency to settle in the dimeric interfaces by direct interaction with the disulfide bridge. These molecules may disturb the dimer formation and may have an inhibition impact on its potential functional role in host immune evasion and virulence pathogenicity. This work provides detailed insights on the sequence and structural variability through computational mutational studies, as well as potent drug candidates with the ability to interrupt the intermolecular disulfide bridge formed between Cys20 and Cys20. Furthermore, the interactions of ORF8 peptides complexed with MHC-1 is studied, and the binding mode reveals that certain ORF8 peptides bind to MHC-1 in a manner similar to other viral peptides. Overall, this study is a narrative of various computational approaches used to provide detailed structural insights into SARS-CoV-2 ORF8 interchain disulfide bond disruptors.

Recent computational drug repositioning strategies against SARS-CoV-2.

Since COVID-19 emerged in 2019, significant levels of suffering and disruption have been caused on a global scale. Although vaccines have become widely used, the virus has shown its potential for evading immunities or acquiring other novel characteristics. Whether current drug treatments are still effective for people infected with Omicron remains unclear. Due to the long development cycles and high expense requirements of de novo drug development, many researchers have turned to consider drug repositioning in the search to find effective treatments for COVID-19. Here, we review such drug repositioning and combination efforts towards providing better handling. For potential drugs under consideration, aspects of both structure and function require attention, with specific categories of sequence, expression, structure, and interaction, the key parameters for investigation. For different data types, we show the corresponding differing drug repositioning methods that have been exploited. As incorporating drug combinations can increase therapeutic efficacy and reduce toxicity, we also review computational strategies to reveal drug combination potential. Taken together, we found that graph theory and neural network were the most used strategy with high potential towards drug repositioning for COVID-19. Integrating different levels of data may further improve the success rate of drug repositioning.

ChikvInt: a Chikungunya virus-host protein-protein interaction database.

Chikungunya is a fast-mutating virus causing Chikungunya virus disease (ChikvD) with a significant load of disability-adjusted life years (DALY) around the world. The outbreak of this virus is significantly higher in the tropical countries. Several experiments have identified crucial viral-host protein-protein interactions (PPIs) between Chikungunya Virus (Chikv) and the human host. However, no standard database that catalogs this PPI information exists. Here we develop a Chikv-Human PPI database, ChikvInt, to facilitate understanding ChikvD disease pathogenesis and the progress of vaccine studies. ChikvInt consists of 109 interactions and is available at www.chikvint.com.

Formal Concept Analysis Applications in Bioinformatics

The bioinformatics discipline seeks to solve problems in biology with computational theories and methods. Formal concept analysis (FCA) is one such theoretical model, based on partial orders. FCA allows the user to examine the structural properties of data based on which subsets of the dataset depend on each other. This article surveys the current literature related to the use of FCA for bioinformatics. The survey begins with a discussion of FCA, its hierarchical advantages, several advanced models of FCA, and lattice management strategies. It then examines how FCA has been used in bioinformatics applications, followed by future prospects of FCA in those areas. The applications addressed include gene data analysis (with next-generation sequencing), biomarkers discovery, protein-protein interaction, disease analysis (including COVID-19, cancer, and others), drug design and development, healthcare informatics, biomedical ontologies, and phylogeny. Some of the most promising prospects of FCA are identifying influential nodes in a network representing protein-protein interactions, determining critical concepts to discover biomarkers, integrating machine learning and deep learning for cancer classification, and pattern matching for next-generation sequencing. [ FROM AUTHOR]

An insight into SARS-CoV-2 membrane protein interaction with spike, envelope, and nucleocapsid proteins.

Intraviral protein-protein interactions are crucial for replication, pathogenicity, and viral assembly. Among these, virus assembly is a critical step as it regulates the arrangements of viral structural proteins and helps in the encapsulation of genomic material. SARS-CoV-2 structural proteins play an essential role in the self-rearrangement, RNA encapsulation, and mature virus particle formation. In SARS-CoV, the membrane protein interacts with the envelope and spike protein in Endoplasmic Reticulum Golgi Intermediate Complex (ERGIC) to form an assembly in the lipid bilayer, followed by membrane-ribonucleoprotein (nucleocapsid) interaction. In this study, we tried to understand the interaction of membrane protein's interaction with envelope, spike, and nucleocapsid proteins using protein-protein docking. Further, simulation studies were performed up to 100 ns to examine the stability of protein-protein complexes of Membrane-Envelope, Membrane-Spike, and Membrane-Nucleocapsid proteins. Prime MM-GBSA showed high binding energy calculations for the simulated structures than the docked complex. The interactions identified in our study will be of great importance, as it provides valuable insight into the protein-protein complex, which could be the potential drug targets for future studies.Communicated by Ramaswamy H. Sarma.

Defective ORF8 dimerization in SARS-CoV-2 delta variant leads to a better adaptive immune response due to abrogation of ORF8-MHC1 interaction.

In India, during the second wave of the COVID-19 pandemic, the breakthrough infections were mainly caused by the SARS-COV-2 delta variant (B.1.617.2). It was reported that, among majority of the infections due to the delta variant, only 9.8% percent cases required hospitalization, whereas only 0.4% fatality was observed. Sudden dropdown in COVID-19 infections cases were observed within a short timeframe, suggesting better host adaptation with evolved delta variant. Downregulation of host immune response against SARS-CoV-2 by ORF8 induced MHC-I degradation has been reported earlier. The Delta variant carried mutations (deletion) at Asp119 and Phe120 amino acids which are critical for ORF8 dimerization. The deletions of amino acids Asp119 and Phe120 in ORF8 of delta variant resulted in structural instability of ORF8 dimer caused by disruption of hydrogen bonds and salt bridges as revealed by structural analysis and MD simulation studies. Further, flexible docking of wild type and mutant ORF8 dimer revealed reduced interaction of mutant ORF8 dimer with MHC-I as compared to wild-type ORF8 dimer with MHC-1, thus implicating its possible role in MHC-I expression and host immune response against SARS-CoV-2. We thus propose that mutant ORF8 of SARS-CoV-2 delta variant may not be hindering the MHC-I expression thereby resulting in a better immune response against the SARS-CoV-2 delta variant, which partly explains the possible reason for sudden drop of SARS-CoV-2 infection rate in the second wave of SARS-CoV-2 predominated by delta variant in India.

HuCoPIA: An Atlas of Human vs. SARS-CoV-2 Interactome and the Comparative Analysis with Other Coronaviridae Family Viruses.

SARS-CoV-2, a novel betacoronavirus strain, has caused a pandemic that has claimed the lives of nearly 6.7M people worldwide. Vaccines and medicines are being developed around the world to reduce the disease spread, fatality rates, and control the new variants. Understanding the protein-protein interaction mechanism of SARS-CoV-2 in humans, and their comparison with the previous SARS-CoV and MERS strains, is crucial for these efforts. These interactions might be used to assess vaccination effectiveness, diagnose exposure, and produce effective biotherapeutics. Here, we present the HuCoPIA database, which contains approximately 100,000 protein-protein interactions between humans and three strains (SARS-CoV-2, SARS-CoV, and MERS) of betacoronavirus. The interactions in the database are divided into common interactions between all three strains and those unique to each strain. It also contains relevant functional annotation information of human proteins. The HuCoPIA database contains SARS-CoV-2 (41,173), SARS-CoV (31,997), and MERS (26,862) interactions, with functional annotation of human proteins like subcellular localization, tissue-expression, KEGG pathways, and Gene ontology information. We believe HuCoPIA will serve as an invaluable resource to diverse experimental biologists, and will help to advance the research in better understanding the mechanism of betacoronaviruses.

Revealing the Molecular Interactions between Human ACE2 and the Receptor Binding Domain of the SARS-CoV-2 Wild-Type, Alpha and Delta Variants.

After a sudden and first spread of the pandemic caused by the novel SARS-CoV-2 (Severe Acute Respiratory Syndrome-Coronavirus 2) wild-type strain, mutants have emerged which have been associated with increased infectivity, inducing surges in the contagions. The first of the so-called variants of concerns, was firstly isolated in the United Kingdom and later renamed Alpha variant. Afterwards, in the middle of 2021, a new variant appeared called Delta. The latter is characterized by the presence of point mutations in the Spike protein of SARS-CoV-2, especially in the Receptor Binding Domain (RBD). When in its active conformation, the RBD can interact with the human receptor Angiotensin-Converting Enzyme 2 (ACE2) to allow the entry of the virions into cells. In this contribution, by using extended all-atom molecular dynamic simulations, complemented with machine learning post-processing, we analyze the changes in the molecular interaction network induced by these different strains in comparison with the wild-type. On one hand, although relevant variations are evidenced, only limited changes in the global stability indicators and in the flexibility profiles have been observed. On the other hand, key differences were obtained by tracking hydrophilic and hydrophobic molecular interactions, concerning both positioning at the ACE2/RBD interface and formation/disruption dynamic behavior.